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Despite the extensive use of in-vitro models for neuroscientific investigations and notwithstanding the 
growing field of network electrophysiology, all studies on cultured cells devoted to elucidate 
neurophysiological mechanisms and computational properties, are based on 2D neuronal networks. These 
networks are usually grown onto specific rigid substrates (also with embedded electrodes) and lack of most 
of the constituents of the in-vivo like environment: cell morphology, cell-to-cell interaction and neuritic 
outgrowth in all directions. Cells in a brain region develop in a 3D space and interact with a complex 
multi-cellular environment and extracellular matrix. Under this perspective, 3D networks coupled to 
micro-transducer arrays, represent a new and powerful in-vitro model capable of better emulating in-vivo 
physiology. In this work, we present a new experimental paradigm constituted by 3D hippocampal networks 
coupled to Micro-Electrode-Arrays (MEAs) and we show how the features of the recorded network 
dynamics differ from the corresponding 2D network model. Further development of the proposed 3D 
in-vitro model by adding embedded functionalized scaffolds might open new prospects for manipulating, 
stimulating and recording the neuronal activity to elucidate neurophysiological mechanisms and to design 
bio-hybrid microsystems. 

Several studies have been devoted to the introduction of in-vitro 3D neuronal systems but the use of such 
experimental models is still limited and, as far as we know, no attempt of functional multisite electro- 
physiological measurements of 3D neuronal networks has been presented in the literature. The potential 
advantages of 3D engineered constructs are evident as they can be used as a more accurate investigational in-vitro 
platform or as the basis for developing living bio-hybrid neuro-electronic microsystems in-vitro or in-vivo'. Thus 
the design and implementation of 3D engineered neuronal networks with embedded sensors and recording- 
stimulating electrodes, would give new opportunities for investigations and applications in the neuroscientific 
domain. On the other hand, the development of such 3D network architectures and the possibility of chronic and 
functional electrophysiological recordings pose new challenges in terms of integration between scaffolds and 
recording-stimulating devices, long-term cell survival, exchange of nutrients, cell coupling with micro-electrodes 
and micro-sensors, etc. Till now, most of the efforts have been devoted to the development of new materials^, 
passive'' '' and active scaffolds^, and new experimental procedures to guarantee the development of 3D cultured 
networks; however, multi-site electrophysiological recordings in such 3D neuronal preparations are stiU lacking. 

Nowadays, the standard and well-accepted experimental model for neurophysiological studies is constituted by 
cultured networks grown onto (stiff) planar substrates also with embedded micro-electrode arrays (MEAs)*". As 
recently pointed out, besides clear advantages related to controllability and observability, such 2D neuronal model 
systems have major limitations, as they might be inherently unable to exhibit characteristics of in-vivo systems. 
For example, soma and growth cones in 2D are unrealistically flattened, and the axons-dendrites outgrowth 
cannot occur in all directions'. Nevertheless, such in-vitro 2D experimental models have been used to investigate 
mechanisms of coding and information transmission', network plasticity and functional connectivity* ', and 
memory'". Although these model systems are somehow widely accepted, one of the raised major issues has often 
been related to the dynamics exhibited by such networks, which are often dominated by bursting activity 
encompassing most of the neurons in the network". Under this perspective the development of a true 3D 
engineered in-vitro neuronal model can certainly be seen as a complementary- alternative and interesting tool 
for neurophysiological investigations. 
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In this work we report a first study about the spontaneous and 
electrically evoked electrophysiological activity in 3D hippocampal 
neuronal cultures coupled to substrate with embedded planar MEAs. 
3D neuronal constructs were developed and implemented by using 
glass microbeads'^ to design engineered in-vitro systems where 
thickness of the preparation, cell density and network connectivity 
are partly controlled to resemble brain tissue while enabling, at the 
same time, optical observation, environmental control and multisite 
electrical recording-stimulation. The culturing method and coupling 
with MEAs is simple, reliable and efficient and it can be easily 
extended to other brain areas (e.g., cortex) and tailored to other types 
of planar MEAs (e.g., high-density micro-electrode array"'") or inte- 
grated with 3D multi-electrodes. 

Results 

We anticipate that the main result is a true 3D neuronal network 
model coupled to planar MEA that shows a wide range of activity 
patterns. These recorded patterns are variable and present striking 
differences from the dynamics of homogeneous-uniform sister cul- 
tures grown onto the planar MEA substrate both during spontaneous 
and electrically stimulated activity. 

3D neuronal network construction. Layers of glass microbeads of a 
selected dimension constitute the basic structure for designing and 
implementing 3D neuronal networks. In order to buUd the 
hippocampal assembly, we used microbeads of 40 |im in diameter 
that self-assembled in a compact hexagonal geometry and that 
provided both a reasonable growth surface for neurites and a large 
enough interstitial space for cell bodies'^. The microbeads were 
sterilized and pre-treated with a mixture of adhesion factors, 
Laminin and Poly-D-Lysine (both 0.05 |ig/ml), at 37°C for 12 
hours. At this stage, about 30' lO"" microbeads were placed onto a 
TransweU® porous membrane where they self-assembled forming a 
uniform layer. For the direct coupling with neurons, we then seeded 
dissociated hippocampal cells (concentration of 700 cells/ |il) on the 
coated monolayer of microbeads (Figure lb, first panel) and cultured 
them in Neurobasal medium (see Methods). The cell seeding 
concentration was chosen in order to obtain a final density of 
about 1,500 cells/mm^ covering the microbeads surface and 
approaching a confluent monolayer. Meanwhile, we plated the 
hippocampal neurons onto the MEA surface, which was previously 
pre-coated with adhesion factors (Laminin and Poly-Lysine both at 
0.05 |ig/ml), in a delimited area obtained by coupling a Poly- 
dunethyl-sUoxane (PDMS) structure to the MEA substrate. The 
confinement (Figure la, first panel) was used to allow cell plating 
preferentially onto the active electrode area of the device and 
subsequently to accommodate for the glass microbeads stacks with 
adherent neurons (Figure lb). The PDMS structure was constituted 
of a simple cylinder (buUt with the molder shown in Figure la, 
second panel), with an external and internal diameter of 22.0 and 
3.0 mm respectively and with a height of 650 |.im. The final 
configuration, composed of MEA device, PDMS structure, 
microbeads and neurons layers, is shown in the third panel of 
Figure la. We wish to underline that the 2D neural network (cell 
density of 1,800-2,000 cells/mm^) directly coupled to the active area 
of the MEA was fundamental to ensure a good communication 
between the substrate embedded micro-electrodes and the 3D 
assembly. 

After 8-10 hours, the suspension of neurons and microbeads was 
moved from the TransweU® membrane and deposited over the 2D 
neuronal network previously plated onto the area defined by the 
PDMS structure (Figure lb, second panel). With respect to the ori- 
ginal density, we estimated a loss of about 20-30% of cells due to 
mechanical manipulations and to cell bodies that remained on the 
surface of the TransweU® membrane. Once the first layer was depos- 
ited, we repeated the same operation to obtain a packed 3D assembly 



(Figure lb, third panel). The resulting 3D structure rearranged itself 
to form a colloidal crystal, which self-assembled in a (theoretically) 
hexagonal structure composed, on average, by 5-8 layers of microbe- 
ads and cells. Once aU the layers were spontaneously assembled, 
neuronal processes grew over the microbead scaffold resulting in a 
high-density physically connected 3D neural network, with an esti- 
mated (on the basis of the actual volume and total number of trans- 
ferred neurons) final cell density of about 80' lO"* cells/mm'', a value 
in the range of the interstitial density computed for glass microbeads 
of 40 |im in diameter (1-10^ mm"^). It is worth noting that this value 
is not far from 92' lO"* cells/mm'', the average neuronal density 
reported for the mouse brain cortex'''. Indeed, by adjusting the 
microbeads dimensions, we could change-control the density of 
the neurons and partly control the average connectivity degree of 
the network'^. 

The final configuration is represented in the cartoon of the last 
panel of Figure lb, where it can be observed that the electrophysio- 
logical activity of the bottom (i.e., read out) layer can be easily 
recorded by the electrodes of the MEA substrate. In addition, the 
developed set-up allowed us to deliver an electrical stimulation both 
from the bottom layer through the MEA electrodes and from the top 
layers by using a conventional tungsten electrode (see Methods). 

To characterize the electrophysiological activity of the 3D network 
(Figure Ic, first panel), and to demonstrate how the 3D structure 
modulates the network dynamics, we specifically compared the 
expressed dynamics with homogeneous-uniform 2D neuronal net- 
works grown over MEAs (Figures Ic, second panel) and with two 
additional controls. The first control (Control 1) was implemented to 
verify whether the electrical activity generated only by 3D networks 
assembled onto the microbeads (i.e., 3D networks coupled to bare 
MEAs) was possibly captured by the substrate MEA (Figure Ic, third 
panel). The second control (Control 2) was constituted of uniform 
2D neural networks plated onto MEAs assembled with around 5-8 
layers (i.e. without neurons) of bare microbeads (Figure Ic, fourth 
panel) to verify whether the electrophysiological activity of 2D net- 
works was affected by the presence of a microbead stack on top of the 
ceU culture. 

Imaging characterization. In order to prove in-vivo like cell 
morphology and demonstrate the reliabUity of implementing a 3D 
network, we gathered structural information on fixed 3D 
hippocampal cultures by indirect immunofluorescence techniques. 
For this purpose we prepared sister cultures on glass coverslips that 
were fixed with PEA 4% (see Methods) together with the 3D-MEA 
cultures at the end of each recording session. We treated the different 
coverslip-grown cultures with a selected panel of primary and 
secondary antibodies in sequence. The use of coverslip-coupled 
(thickness 0.13-0.14 mm) 3D cultures was necessary to overcome 
the limitation introduced by the relatively high thickness (~1 mm) 
of the MEA glass substrate that precluded the possibility of using 
high numerical aperture objectives with an inverted confocal 
microscope. The configuration of the 3D-coverslip cultures was 
identical to the final structure of the hippocampal 3D-MEAs. 
Hippocampal neurons were plated by the same procedures (see 
Figure lb and Methods) and for the same period. We also used an 
up-right confocal microscope to complete our tridimensional optical 
analysis. The objective was immersed inside the MEA reservoir 
where the multilayered structure of microbeads was assembled. In 
this scenario, we used long-working distance objectives, to observe 
the morphology of the neural network within the deepest layers 
with a z-stack exceeding 400 |j,m. The pictures in Figure 2 show 
the different configurations investigated: from 2D cultured hip- 
pocampal neurons to 3D networks coupled to MEAs (4* week 
in-vitro). Particularly, Figure 2a shows a close-up of a few ceUs 
cultured on a planar substrate, and Figure 2b shows a small cluster 
of neurons with embedded microbeads in which a GABA and NeuN 
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3D network on MEA 2D network on MEA Control 1 Control 2 

3D network on bare MEA 2D network with 

bare micro-beads 



Figure 1 | Construction of 3D nevu-al networks, (a) . From left to right: PDMS structure allowing to confine microbeads and neurons on the recording site 
area; molder used to build the confinement structure; cartoon that illustrates the final configuration (multi-layers of microbeads and neurons confined by 
a PDMS structure onto the active area of the Micro-Electrode Array), (b). Main steps for building a 3D neural network. Microbeads were placed onto a 
porous Transwell® membrane where they self-assembled in a hexagonal geometrical structure; dissociated hippocampal cells were plated on such coated 
microbeads. To obtain a 3D structure, the suspension of neurons and microbeads was then moved from the membrane to the MEA surface several times 
(details can be found in the Results Sec). The last sketch depicts the recording/stimulation configuration: the electrophysiological activity of the 3D 
network is recorded from the substrate MEA electrodes (bottom layer); the network is stimulated by using both MEA and tungsten electrodes (bottom 
and top layers), (c). Experimental configurations: 3D neural network; 2D neural network; Control 1 i.e., network constructed among beads plated on 
MEA; Control 2 i.e., 2D network plated onto MEA with bare beads on top. 



positive neuron (top left) and four NeuN positive (putative 
excitatory) neurons can be observed. It is worth noting the round 
shape of the cell bodies with respect to the flattened shape of the 
surface plated neurons. Panels c-f of Figure 2 show a sequence of the 
microbead stack and neurons (also visible as Supplementary - Movie 
1) on a glass coverslip. On the orthogonal projection of the y-axis it 
can be observed some cell bodies lie on the glass surface while others 
remain entrapped among the microbeads in the interstitial spaces 
(see also Supplementary - Movie 2). Figure 2g shows extensive 
encapsulation by neuritic processes (MAP-2 positive) around the 
microbead scaffold indicating a high-density of potential con- 
nectivity among neurons. Finally, Figure 2h and 2i show the floor 
of the 3D networks (with the substrate embedded microelectrode) 
and a nice occupancy of the interstitial space and extensive 
arborizations connecting the neurons. In general, we can underline 
a main difference in the development of the axo-dendritic processes. 
The growth of 3D network shows patterns of tangled ramifications 
both across and along different layers (5-8 layers, see Supplementary 

- Movie 3), and around the single microbead. This new network 
configuration gains a wide range of extension (see Supplementary 

- Movie 4) compared to the 2D model and, although constrained by 
the presence of the microbeads ordered reticulum, each neuron, in 
principle, can have a site density of connecting points that is inversely 



dependent on the beads radius cubed'^. Considering the close 
relationship between the spontaneous electrophysiological activity 
and the level of synaptic expression"", we can assume that the large 
number of synaptic puncta visible in Figure 2c-f (see also 
Supplementary - Movie 1) is an indication of a possible high 
connectivity degree wdthin the new network model. From these 3D 
neural networks we could estimate a synaptic density of about 
2.6' 10^ synapses/mm' (see Supplementary - Fig. SI, a-c), value 
that is an order of magnitude less than what is reported for the 
mouse cortex'^. By considering the average neuronal density of the 
3D model (80* 10' cells/mm') and the estimated synaptic density, we 
obtained about 600 synaptic connections per neurons. As a reference, 
we also estimated the average number of synapses per neuron of the 
considered 2D neuronal cultures (see Supplementary - Fig. SI d-e) 
obtaining also in this case an average density of —600 synapses/ 
neuron". 

Spontaneous activity. Spontaneous activities of 2D (n = 28) and 3D 
(n = 39) neural networks were recorded for 30 min after 3-4 weeks 
in culture. Comparisons with control experimental cultures (3D 
cultures on bare MEA - Control 1, and 2D cultures coupled to 
bare 3D microbeads structures - Control 2) were also carried out, 
to verify possible unwanted effects not attributable to the 
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Figure 2 | Imaging characterization, (a). Mature hippocampal neurons (4th week in-vitro), cukured on a planar substrate (i.e., 2D standard model, 
positive for Map-2 andDapi nuclear dye). Scale bar: 20 |j.m. (b). Small cluster of neurons and microbeads closely linked together by a GABA interneuron 
(green/red) and four Nuclear Neuronal protein NeuN positive neurons (red). Scale bar: 20 |j.m. (c-f). z-stack sequence of 3D cultured hippocampal 
neurons fixed at DIV (Days In Vitro) 28 and immunolabeled for the pre-synaptic marker Synapsin (green) and NeuN (red). Scale bar: 10 |j.m. (g). 3D 
rendering of 26 1 .4 \rca z-stack of the hippocampal network on MEA device showing a high density of neuronal processes, exposed to dendritic marker 
MAP-2. Scale bar: 40 ^m. (h). 3D culture on MEAs, labeled with MAP-2 and NeuN. Scale bar: 40 |j.m. (i). The first layer of microbeads directly coupled to 
the top of a single micro-electrode. Scale bar: 20 |im. Note: secondary antibodies for all images are Alexa Fluor 488 or Alexa Fluor 546 conjugated with 
Goat anti mouse or Goat anti rabbit. 



development of a true 3D network. All cultures were initially 
investigated with both excitation and inhibition active. The role of 
inhibition in network dynamics is treated later. 

First of all, the spontaneous activity of 3D hippocampal neur- 
onal networks was successfully recorded in all the reported experi- 
ments by means of planar TiN 60 channels MEAs (see Methods). 
The electrophysiological activity of the 3D networks was then 
analyzed and compared with the uniform 2D hippocampal cul- 
tures, and with the two additional control conditions (Figure Ic). 
As reported in the literature^'"''", 2D uniform hippocampal and 
cortical neuronal networks show a typical electrophysiological 
activity in the mature phase of their development (i.e., after 3-4 
weeks in-vitro) characterized by oscillating dynamics and syn- 
chronous bursts involving most of the active channels in the 
MEA. Figure 3a shows 5 min of the spontaneous activity of a 
representative experiment performed on a 2D hippocampal culture 
during the fourth week in-vitro. Also in our recordings, we found a 
quasi-synchronous activity and a network dynamics composed 
mainly of network bursts (NB) (i.e., synchronized events that 
involve most of the recording channels; see Methods). Looking 
closer at the raw data of a single channel, one can see that these 



bursts are characterized by an amplitude of about 100 nV and a 
duration of about 300 ms (see also Supplementary - Movie 5). 
Similarly, Figure 3b shows a raster plot of the spontaneous activity 
of a representative experiment of a sister culture (i.e., same batch) 
forming a 3D network at the same age of its development (i.e., 
fourth week). It can be noticed that the signature of the network 
dynamics presents a wider repertoire of activities with less global 
synchrony and more pronounced spiking at a single channel- 
neuron level. From a preliminary qualitative observation, we can 
see periods in which the network holds a more synchronous activ- 
ity (red box) similar to what observed in the 2D networks. In this 
case bursts are characterized by a similar amplitude (i.e., about 
100 nV), but by a longer duration (up to 1 s, see raw data of a 
single channel displayed in figure 3b top; see also Supplementary - 
Movie 5). As anticipated, the 3D neural network exhibits also a 
significant 'random spiking' and non-synchronous bursting activity 
(blue box). To highlight these two distinct 'modes' of firing, we 
computed (cf. boxes of Figure 3b) the Instantaneous Firing Rate 
(IFR) profiles (see Methods). The same plot for the 2D network 
(Figure 3a), relative to the green region, shows a very high synchron- 
ous activity and the absence of significant random spiking activity. 
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Figure 3 | Spontaneous activity characterization. Electrophysiological activity of a representative neuronal network, (a). 2D and (b). 3D. From top to 
bottom: example of 10 s of raw data recorded from one microelectrode; raster plots showing 300 s of spontaneous activity. For the highlighted areas, 
the Instantaneous Firing Rate (IFR) profile is shown. For the 3D network, the blue and red boxes highlight random spiking and network bursting activity 
regions, respectively. For the 2D network, only network bursting regions can be detected (green box), (c). Percentage of random spikes, (P-value = 10"*, 
Mann-Whitney U-test). (d). Mean Bursting Rate (P-value = 10"*, Mann-Whitney U-test). (e). Burst Duration, (f). Mean Network Bursting Rate 
(P-value = 9- lO"*, Mann-Whitney U-test). (g). Network Burst Duration (P-value = 9- 10"*, Mann-Whitney U-test). (h). Inter Burst Interval 
distribution with a bin size of 0.5 ms. Parameters were evaluated over the entire dataset. Asterisks above the plots indicate statistically significant 
differences. 



The previous observations were further globally quantified by 
evaluating the percentage of random spikes (i.e., the fraction of spikes 
outside bursts) in all the performed experiments (Figure 3c). 2D 
networks exhibit a strong synchronous bursting activity with a low 



percentage of random spikes (15.9 ± 2% (mean ± standard error of 
the mean)), differently the 3D networks show a significant higher 
level of random spiking (33.8 ± 3%; P-value = 10"*, Mann- Whitney 
U-test). In order to quantify the bursting activity at the single channel 
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level, the frequency and the duration of the bursts were evaluated in 
all the performed experiments (Figure 3d-e respectively). 2D net- 
works presented a much higher level of bursting rate (19.7 ± 
1.1 burst/min) than 3D neural networks (5.2 ± 0.4 burst/min, P- 
value = 10"'', Mann- Whitney U-test), while burst durations are not 
significantly different between the two. Even if network bursting 
remains a common signature for both preparations, the shift from 
2D to 3D did also affect the bursting behavior as the main effects are 
observed in both frequency and duration. Then, we computed the 
network mean bursting rate and duration evaluated over the entire 
dataset (Figure 3f-g respectively). In 2D neural network the network 
mean bursting rate is high (22.3 ± 3.8 NB/min) with a nearly con- 
stant duration (540 ± 5 ms); differently, the 3D networks exhibit few 
network bursts (3.0 ± 0.6 NB/min, P-value = 9-10"*, Mann- 
Whitney U-test) with a high duration (919 ± 65 ms, P-value = 
9" 10"'', Mann- Whitney U-test). The high value of the standard error 
for the 3D networks indicates a much higher variability of bursting 
activity. Further, we evaluated the Inter Burst Interval (IBI) distribu- 
tions for both 2D and 3D neural networks (Figure 3h). For the 
uniform 2D networks, about 90% of the burst intervals were distrib- 
uted within the first 2.5 s and only 1% of the burst showed IBIs 
greater than 5 s. For the 3D networks, IBIs are much shorter and 
about 40% of the burst intervals are in the first bin (0.5 s). The 
remaining distribution is much broader and covers intervals up to 
10-15 s. We also challenged the presented 3D model with the addi- 
tion of bicuculline (BIC, 25 |,iM) to test whether, by blocking GABA 
transmission, we had a different effect on 3D networks with respect 
to uniform 2D cultures. We found that, starting from different spon- 
taneous dynamics, the administration of BIC steered both 2D and 3D 
neuronal networks towards the same more synchronized dynamic 
state (see Supplementary - Figure S2) with a similar percentage of 
random spike and NB duration (Figure S2d and h). Interestingly, also 
the IBI distributions become similar after the addition of BIC. 

We also quantified the level of synchronization and correlation of 
the network activity by means of the cross-correlation (CC) function 
(see Methods) computed for 2D and 3D neural networks with a time 
window of 2 s and a bin size of 5 ms. For 3D neural networks, CC was 
also specifically evaluated within the random spike activity and 
within the burst activity regions. CC function for a 2D network 
(Figure 4a) underlines a high level of correlation among a represent- 
ative channel and all the other recording channels. The Cpeak values 
(i.e., maximum values of the CC function) for all the active channels 
are close to one (0.85 ± 0.03). On the contrary, for a 3D network we 
can see a low level of correlation among a representative channel and 
all the other active channels in the random spike activity region 
(Figure 4b). The mean Cpeak value is low (0.42 ± 0.02) and the peaks 
values are shifted with respect to the central bin suggesting possible 
causal relationships mediated by the neuronal populations above the 
read-out layer. Within the burst activity region, the level of correla- 
tion among channels is high (mean Cpeak value of 0.70 ± 0.02), and 
the majority of peaks is centered in the central bin (Figure 4c). We 
further characterized the level of synchronization of the 2D and 3D 
network activity by computing the values of the CC function in the 
central time bin (i.e., C(0)). 2D neural networks exhibit a high level of 
synchronization as the normalized C(0) values are close to one for the 
majority of the recording channels (Figure 4d). The level of syn- 
chronization of the 3D neural network activity is high during the 
period of bursting (C(0) close to one as for the 2D networks - 
Figure 4f) while it is very low (Figure 4e) within the random spike 
activity region. In order to better identify the synchronized channels 
within this region, we compared the C(0) values with a threshold 
related to the average C(0) values (see Methods). Only five channels 
presented a value of C(0) higher than the threshold (these channels 
are highlighted in pink in the raster of Figure 3b) and they represent a 
group of spatially confined electrodes that could suggest the presence 
of a synchronized sub-population. Putative assembly of neurons can 



be individuated by higher C(0) values shown in the 8X8 color map 
layout which represents also the physical position of the MEA chan- 
nels (Figure 4g-i). In the 2D neural network several channels present 
a maximum C(0) value; in the 3D network (within the random spike 
activity regions; 3D-rs) there is only one channel with a maximum 
C(0) value (autocorrelation) and there are few channels with a C(0) 
value greater than 0.5. In the 3D burst activity regions (3D-b), beside 
the autocorrelation, almost all the other channels exhibit a C(0) value 
greater than 0.7. The previous observations were globally quantified 
by evaluating the C(0) values for all the performed experiments 
(Figure 4j). 2D neural networks present a mean C(0) value close to 
one (0.81 ± 0.02), suggesting a very high synchronization among all 
the channels. The level of synchronization in 3D neural networks is 
high within the burst activity region (mean C(0) value of 0.7 1 ± 0.04) 
and low within the random spike activity region (mean C(0) value of 
0.46 ± 0.02). Data differences are significant (P-value < 10"*, Mann- 
Whitney U-test). To prove such activity variability was supported by 
the development of a true 3D architecture, we compared the exhib- 
ited 3D dynamics with the activity of the Controls presented in 
Figure Ic. As it might be expected, in the case of Control 1 (n = 
5), almost no activity was recorded (see Supplementary - Fig. S3a). 
This behavior is due to the assembly procedure by which neurons are 
only seeded within the microbead structure, thus resulting in a bad 
coupling of the cell bodies to the MEA substrate (see Figure Ic). It 
should be underlined that only with a good contact (i.e. effective 
adhesion) between neurons and microelectrodes we can reliably 
record electrophysiological signals''^. In the case of Control 2 (n = 
5), the electrophysiological activity was, as expected, very similar to 
the 2D neural network (see Supplementary - Fig. S3b). Moreover, we 
compared the activity of Control 2 with the activity of uniform 2D 
neural networks, by evaluating the percentage of random spikes, the 
NB duration and the C(0) values for all the performed experiments 
(see Supplementary - Fig. S3c, S3d, S3e). These results clearly indi- 
cate that the presence of a 3D network organization introduces a 
great variability in the electrophysiological activity recorded at the 
bottom layer, by inducing partly synchronized sub-populations and 
breaking down the uniform and stereotyped behavior of uniform 
high density 2D neuronal networks^". 

To demonstrate whether the layers of the 3D structures were 
effectively connected and capable of inducing significant differences 
in the recorded dynamics, we studied and characterized the propaga- 
tion of the electrophysiological signals in the read-out layer (i.e., the 
bottom layer of the neurons coupled to the planar MEA). In the 2D 
neural networks, the signal propagation is indubitably continuous 
(see Supplementary - Movie 5). Bursts start from few sources and 
propagate rapidly and in a continuous way invading the whole net- 
work, as also observed in calcium imaging recordings^'. In the 3D 
networks, the signal that starts from one site looks as if it propagates 
in the network in a discontinuous way (see Supplementary - Movie 
5). The observed variability measured at the read-out layer can be 
explained by considering a propagation also in the z-direction. To 
confirm the impression of activity propagation also in the upper layer 
of the 3D networks, we evaluated the ignition sites from which the 
NBs started. We defined such ignition sites as the sites where bursts 
initiate in at least 5% of all the NBs^^. Moreover, we accounted for the 
latency of the subsequent spikes in all the recording electrodes (i.e., 
the time between consecutive spikes in another electrode involved 
into the NB starting from the same source). Figure 5a shows the 
occurrence of single sources for 30 min of activity of a representative 
2D network. It is worth noting that only three major sources were 
individuated. Furthermore, such signal sources were stable in time as 
we evaluated the occurrences of these three electrodes by dividing the 
recorded interval in phases of 5 min each. The three sources were 
present in all the four considered experimental phases (Figure 5b). 
Figure 5c shows the map of the mean latency (interpolated propaga- 
tion fronts, isochrones) between a chosen electrode (source) and the 
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spikes which follow in all the electrodes involved in the NBs initiated 
by that electrode^"". Following the same approach, we analyzed the 
signal propagation during NBs in a representative 3D network. We 
found a large number of sources with low occurrences (Figure 5d) 
and by evaluating the signal sources in consecutive 5 min experi- 
mental phases, we found a great variability (Figure 5e) as different 
sources are seemingly detected in different phases. The isochrones 
map (Figure 5f) displays a discontinuous propagation of the signals 
with long time delay. The previous illustrative results were further 
globally quantified by evaluating the number of NBs and the relative 
detected sources in six different experiments (Figure 5g). The num- 
ber of NB exhibited by the 2D networks was higher (thick gray bars) 
compared to the number of NB in 3D networks. Indeed, we found 
more ignition sites (thin black bars), in the 3D networks than in the 
2D ones. Summarizing, 2D networks exhibited a lot of NB with few 
sources; 3D networks few NBs with more seeming sources (Figure 5g, 
inset). We also looked at possible correlations between latency and 
distance. For 2D networks, we found a linear relationship between 
distance and latency for all the analyzed experiments (R^ = 0.67; 
Figure 5h, black line). The estimated speed of propagation was 42 
± 3.5 mm/s, value comparable with previous works^''^'*'^"'. From the 
3D networks (Figure 5i), we could not extract a good correlation 
between latency and distance (R^ = 0.07; black dashed line) meaning 
that only random latencies (up to 100 ms), with respect to the dis- 
tance, were observed in the read-out layers of neurons coupled to the 
planar MEA. 

The discontinuity in the propagation and the seeming variability 
of the signal sources, are all indirect proofs that signals propagate in 
the upper layers and support the hypothesis of the presence of signal 
sources in the upper layers. 

Stimulus-evoked activity. To further support the previously pre- 
sented results and to ultimately demonstrate whether the layers were 
functionally connected (i.e., synaptic pathways cross the assembly in all 
the directions of the 3D space), we performed experiments by 
delivering electrical stimulation with two different experimental 
protocols-settings. We stimulated from the bottom layer 3D (n = 
15) and 2D (n = 7) cultures by using one of the electrodes of the 
substrate MEA; we performed top-layer stimulation on 3D cultures (n 
= 3) by means of a bipolar tungsten electrode. The experimental 
configuration is illustrated in the last panel of Figure lb. 

Figure 6 shows the Peri Stimulus Time Histograms (PSTHs) over a 
4X4 grid (a subset of the channels of the 8 X 8 MEA layout) with the 
different evoked activities that the network exhibited in response to 
the applied stimulation. To examine the network response charac- 
teristics, we evaluated the time in which the 60% of the PSTH area 
occurred (tgo). The vertical red bars in Figure 6 indicate such tgo 
value. In the 2D network, the neural response to the stimulation 
delivered by the bottom layer (2D-Bottom) was, as expected^"", sim- 
ilar for all the responding channels (Figure 6a). The highest response 
peak (i.e., the maximum of the evoked response) was detected imme- 
diately after the stimulation; then the number of evoked spikes exhib- 
ited a decreasing trend that goes to zero in about 1 s. The 60% of the 
evoked spikes (see red bar, tgo) was confined to the first hundreds of 
milliseconds. In the case of the considered 3D networks stimulated 
by the MEA (Figure 6b, 3D-Bottom), the PSTH profiles presented a 
more pronounced variability. Some electrodes displayed a large res- 
ponse immediately after the stimulation. In other cases, the response 
peak was up to 1 s after the stimulation. In the highlighted case, the 
single channel PSTH profile indicated the 60% of response was 
within the first 700 ms. These results indirectly prove by themselves 
that additional synaptic pathways were involved and elicited when 
stimulated from the bottom layer (MEA substrate) and that possible 
reverberating paths going up and down into the 3D microbead layers 
contributed to the observed delayed responses. The last analyzed case 
concerns the 3D network response to a stimulation delivered by the 



top layers (Figure 6c, 3D-Top). As it might be expected, in this 
configuration not many channels (less than 50%) detected a neural 
response, suggesting that a reduced number of synaptic pathways 
reached the bottom layer. This could be an indication that the con- 
nectivity along the z- direction was less effective or that the synaptic 
density decreased from bottom to top. However, for the responsive 
channels, the evoked activity was quite different from the previous 
ones. The tgo (i.e., red bars) spanned over a wide time window and 
electrodes exhibited a tgo value exceeding 1 s. Finally, no response 
was recorded in the Control 2 condition from top stimulation, indi- 
cating there was no wide field effect of the stimulation. 

In the illustrative case of Figure 7, for the 2D network stimulated 
from the bottom layer, the responses were close to the stimulation 
time (Figure 7a) and, for the majority of the channels, the 60% of the 
PSTH area (i.e., t^o) was detected within the first 250 ms (Figure 7d) 
with all the electrode responding the same way. For the 3D neural 
networks stimulated from the bottom layer, there was a high vari- 
ability of the tgo and such parameter shows an almost uniform dis- 
tribution within an interval of 100-1,200 ms (Figure 7e). Finally, in 
the case of 3D networks stimulated from the top layer, the responses 
to stimulation were very much delayed (Figure 7c) as the recorded 
evoked responses (related to neurons in the read-out layer) needed to 
be transmitted to such bottom layer with the possible involvement of 
several neurons and reverberating circuits. The t^o time values dis- 
tribute in an interval ranging from 500 ms to 1,400 ms (Figure 7f). 
The cut-off value observed around 500 ms might indicate that only 
very few poly-synaptic pathways or direct functional connections 
crossed this culture in the z-direction as the 60% of the total res- 
ponses reached the MEA substrate after such a long time delay. The 
electrophysiological specificity of the analyzed 3D networks was fur- 
ther demonstrated by the analysis and comparison, in all the per- 
formed experiments, of the tgo and tgg (time delay to have the 90% of 
the response to the stimulation) (Figure 7g). Time values are reported 
for 2D and 3D networks stimulated from the bottom layer (P-value 
= 9*10"'', Mann- Whitney U-test) and for 3D networks stimulated 
from the top layer (P-value = 2' 10"'', Mann- Whitney U-test). 

2D networks displayed a reliable response time (also qualitatively 
found in previous works'""'^') in which the 60% of the PSTH area 
occurred within the first 200 ms. The values detected for the 3D 
networks were higher and more variable. In the case of the bottom 
layer stimulation, the tgo value was 680 ms ± 28 ms, and in the case 
of the top layer stimulation, the tgo shifted to 1,050 ms ± 38 ms. The 
time in which we found the 90% of the total evoked responses to 
stimulation (tgo) was obviously higher than tgo, but it followed the 
same trend (Figure 7g). 

Moreover, we evaluated the percentage of the PSTH area gener- 
ated after electrical stimulation in different time windows (Figure 7h) 
to quantify the different response modes of the measured networks. 
The response to stimulation of 2D networks was very fast: about 90% 
of the PSTH area was reached within the first 500 ms. 3D neural 
networks stimulated from the bottom layer exhibited a similar res- 
ponse to stimulation within all the time windows: it can be observes a 
fast response (35% of PSTH area reached within the first 100 ms) and 
also a slow response (about 30% of the response to stimulation was 
detected between 1 s and 2 s). In the case of 3D networks stimulated 
from the top layer, the fast response was not detected (only about 
10% of PSTH area was within the first 300 ms) suggesting that no 
detectable direct pathway was present. The responses generated from 
this type of stimulation were slow: about 90% of the PSTH area was 
reached between 500 ms and 2 s. 

As a final control, to prove that evoked responses recorded in the 
read-out layer when stimulating from the top, were due to synaptic 
transmission, we performed specific experiments by blocking syn- 
aptic transmission. We administered a cocktail of bicucuUine (BIC, 
25 |.iM), 2-amino-5-phosphonovalerate (APV, 50 |.iM), 6-cyano-7- 
nitroquinoxaline-2,3-dione (CNQX, 50 |xM)^'' after a first phase of 
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repeated stimulation (control condition) of the analyzed 3D network. 
Before the synaptic blockade, the top layer stimulation indeed gen- 
erated a response of the network (see Supplementary - Fig. S4). When 
we administered the synaptic blockade to the culture, we did not find 
any evoked response (see Supplementary - Fig. S4b-c). This result 
indicates that the response detected at the bottom layer was due to 
the synaptic transmission among different neurons that eventually 
found a target neuron coupled to the planar MEA in the read-out 
layer. Thus, the evoked response was due to synaptic reverberation 
and amplification and not to purely neuritic propagation meaning 



that direct connections from top layer to the surface of the MEA are 
unlikely. That is also somehow intrinsically due to the experimental 
procedure for the 3D network construction in which the bottom 
layer is first coupled to the planar MEA and then assembled with 
the upper layers. 

Discussion 

We presented a new experimental in-vitro platform constituted of 
3D engineered neuronal cultures coupled to MEAs for network 
electrophysiology. 
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The use of microbead scaffolds was specifically optimized and 
tailored to be integrated with planar MEA in order to investigate 
the functional properties (i.e. electrical activity) of 3D hippocampal 
networks that were never presented before. We carefully compared 



the spontaneous and evoked electrophysiological activity of 3D 
cultures with the traditional network dynamics from homogen- 
eous-uniform 2D cultures. First, in the developed 3D engineered 
hippocampal assemblies, cell morphology, connectivity (see also 
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Figure 2) and, potentially, extracellular matrix, are much closer to the 
in-vivo situation' *'^". Second, we observed striking differences with 
respect to uniform 2D networks, and we quantified evidences of a 
wide dynamic repertoire of activity patterns possibly recapitulating 
various activity dynamics of in-vivo brain regions^''''". Quasi- 
synchronous network bursting activity of high density 2D neuronal 
networks'^'"*"''^ was maintained in 3D cultures but global frequency 
was decreased (<1 Hz) and global synchrony was lost (spatial 
segregation of bursting - see Figure 3). Even if the recorded dynamics 



were very different, the average connectivity (i.e., 600 connections 
per neuron) was found to be the same in both 2D and 3D networks 
and the proportion between GABAergic and glutammatergic popu- 
lations was likely to remain similar in both 2D and 3D models (with 
about a 1 : 5 ratio, see also Supplementary - Fig. S5 a-b)^'*'^''"'''. In the 
3D architecture, we can speculate that, with the enhanced dendritic 
development, the resulting network was more complex, and that 
cell morphology might also contribute to such increased complexity. 
We can hypothesize that the different GABAergic interneuron 
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populations could gain a considerable increase of the functional 
surface (volume) compared to those grown in the 2D cultures. 
These possible wider and longer interactions with other neurons 
could contribute to desynchronize and temporally differentiate the 
network activity, hence producing bursting activities confined to 
sub-populations and global asynchronous activity patterns^^. The 
addition of BIC induced in the 3D network model a more synchro- 
nized state (similar to what happen in uniform 2D networks), indi- 
cating a possible more synchronized activity in the upper layers, and 
thus determining a reduced impact of the electrophysiological signals 
coming from above (see Supplementary - Fig. S2). 

Along this direction, a systematic study of the spontaneous activity 
of 3D networks during development with the specific labeling of 
interneurons, could contribute to the understanding of the signifi- 
cance of the balance between excitation and inhibition with specific 
reference to the role of GABAergic populations^^^"^^, distribution of 
hubs^*^ and network synchronizations^'^. 

Although such variable spontaneous activity together with the 
analysis of signal propagation is a convincing hallmark of the inher- 
ent contribution of the 3D structure to the network dynamics, the 
ultimate evidence that we are playing with actual 3D functionally 
interconnected neuronal networks is given by the results obtained 
upon the application of an external electrical stimulation. 

2D interconnected networks stimulated by different electrodes, 
showed a relatively fast and synchronized response depending on 
the stimulating site. The site of stimulation causes an almost uniform 
response (for the high-connectivity of the network) in all the active 
sites (see Supplementary - Fig. S6a-b). The response time is rela- 
tively fast and confined below the first 250 ms (Figure 7d). 3D net- 
works stimulated from the bottom layer exhibit various types of 
synchronized response (see also Supplementary - Fig. S6c-d) often 
generating a uniformly distributed response delay (Figure 7e). In 
general, the evoked activity can be thought as the combination of 
two types of response. The first type is generated by subsets of neu- 
rons in the read-out layer that are connected with subgroups of 
neurons in the same layer. The second type of response is generated 
by subsets of neurons in the read-out layer connected with sub- 
networks of neurons in the upper layers. The first type of response 
is fast; the second type is slow. Because of the different described 
types of activated circuits, the network response is not globally syn- 
chronized. When stimulating 3D networks from the top layer, the 
stimulation needs to be transmitted, through different neurons of the 
various layers, to neurons of the read-out layer coupled to the MEA 
generating time delays greater than 500 ms (cf. Figure 7f). These long 
delays could be explained by a less effective connectivity in the z- 
direction and by a possible decrease in the synaptic density from the 
bottom to the top layer. The administration of synaptic blockers to 
the culture impairs the recording of evoked responses in the bottom 
layer, thus demonstrating that activity propagates through synaptic 
reverberation and amplification in the different layers. The enhanced 
variability in the time delays of the evoked responses (cf Figure 7h) is 
an additional proof we are facing heterogeneous functionally inter- 
connected 3D cultured networks. 

In summary, the developed experimental paradigm can constitute 
the basis for a next class of experimental models to study neurophysi- 
ology in in-vitro systems, for investigating the computational prop- 
erties of neuronal networks, and for developing new bio-hybrid 
microsystems. 

Methods 

Cell culture. Hippocampal neurons were dissociated from E18 Sprague Dawley rats 
(Charles River Laboratories, Milano). The procedure was approved by the European 
Animal Care Legislation and by the guidelines of the University of Genova. The day 
before the plating, MEAs {Multi Channel System MCS GmbH) were sterilized at 
120 'C in the oven, and the microbeads were sterilized with ethanol at 70% for 1-2 
hours. The sterilized material was therefore exposed to the coating treatment with 
adhesion protein Laminin (L-2020 Sigma) and Poli-Lysin (P-6407 Sigma) at 0.05 [ig/ 
ml overnight in incubator at 37 C. The adhesion factors were removed and the treated 



glass surfaces washed with sterile water. MEAs were left to dry inside the laminar 
hood. Lastly, the PDMS mask (Corning Sigma) was put on the electrodes area. From 
each embryo, hippocampi were removed and placed into ice cold Hank's balanced salt 
solution. The tissue was then dissociated in 0.125% of Trypsin/Hank's solution 
containing 0.05% of DNAse (D-5025 Sigma- Aldrich) for 15-18 min at 37 C. The 
supernatant solution was removed and the enzymatic digestion was stopped by 
adding 10% fetal bovine serum (EES) in Neurobasal medium for 5 min. Medium with 
EES was removed and replaced with Neurobasal medium supplemented with B27, 1 % 
Glutamax, gentamicin 10 [ig/ml (Gibco Invitrogen). Some cells were then plated 
directly onto the area of MEA electrodes, defined by the PDMS mask, while some 
others were distributed onto the microbead monolayer inside the Transwell® 
supports (Costar Sigma) to create the 3D assembly. The obtained 2D and 3D cultures 
were maintained in a humidified CO2 atmosphere at 37' C for 4-5 weeks. Half of the 
media was replaced once a week. No antimitotic drug was added to prevent glia 
proliferation, since glial cells are essential elements for a healthy development of 
neuronal networks. 

Glass microbeads. Glass microbeads (distributed by DistrUab-Duke Scientific) were 
used as scaffold to create 3D neuronal networks. The cell density and the thickness of 
the reconstructed tissue can be tuned depending on the diameters of the sphere and 
the number of transferred layers. Glass microbeads with diameter of 40 |j.m (certified 
mean diameter of 42.3 ±1.1 |,im) were used and 5 to 8 layers were stacked. 

Immunocytochemistry. Hippocampal cultures were fixed in 4% paraformaldehyde 
in Phosphate Buffer solution (PBS), pH 7.4 for 30 min at room temperature. 
PermeabUization was achieved with PBS containing 0.5% Triton-XlOO for 15 min at 
room temperature and non-specific binding of antibodies was blocked with an 
incubation of 45 min in a blocking buffer solution consisted of PBS, 3% BSA (bovine 
serum albumin Sigma) and 0.5% FBS. Cultures were incubated with primary antibody 
diluted in PBS Blocking buffer for 2 hours at room temperature or incubated at 4°C 
overnight in a humidified atmosphere. Cultures were rinsed three times with PBS and 
finally exposed to the secondary antibodies: Alexa Eluor 488 or Alexa Eluor 546 Goat 
anti mouse or Goat anti rabbit, diluted respectively 1 : 700 and 1 : 1,000 (Invitrogen 
Life Technologies S. Donato Milanese). The antibodies used in the images are in the 
following order: MAP-2 1 : 500 (Chemicon Millipore), NeuN 1 : 400 (Chemicon 
Millipore), GABA 1 : 700 (Sigma), Synapsin 1 : 200 (Chemicon Millipore), Dapi 
300 nM (Sigma). To verify the presence of glial cells in the culture, we fixed and 
exposed to the marker GEAP (monoclonal or polyclonal antibodies Sigma) both 2D 
controls (on glass and on MEA) and 3D samples assembled on coverslips or onto 
MEA. The evaluation of the samples was performed by confocal-laser scanning 
microscopy, widefield epifluorescence microscope Olympus IX- 70 inverted (Fig. 2a) 
and Olympus BX-51 upright microscope (Eig. 2b). Image acquisition with Olympus 
was done with a Hamamatsu Orca ER II digital cooled CCD camera driven by Image 
ProPlus software (Media Cybernetic). 

Optical setup. Confocal imaging was acquired on two different microscopes: a Leica 
TCS SP5 AOBS Tandem DMI6000 inverted microscope using a HCX PL APO 63 X 
N.A. 1.40-0.60 Oil objective (Leica Microsystems, Mannheim, Germany) and a Nikon 
Al Confocal inverted microscope equipped with CElPlan Apocr XlOx. W.D.I mm 
A.N. 0.75 (Nikon Instruments Europe B.V, Amsterdam). 

2PE imaging was performed using a fs-pulsed Ti: sapphire laser Chameleon Ultra 
(Coherent, Santa Clara, CA, USA) coupled with Leica TCS SP5 AOBS Tandem 
DM6000 upright microscope equipped with HCX APO 20X N.A. 0.50 WI objective 
and HCX APO L 40X N.A. 0.80 WI objective (Leica Microsystems, Mannheim, 
Germany). Data were analyzed by means of the LAS AE V.2.6 software (Leica 
Microsystems, Mannheim, Germany), Nis elements V4.0 (Nikon Instruments 
Europe B.V, Amsterdam), ImageJ (National Institutes of Health, USA) and SVI 
Huygens Professional for volume rendering. 

Data analysis. Data analysis was performed off-line by using a custom software 
package named SPYCODE^ developed in MATLAB (The Mathworks, Natick, MA, 
USA). Spike detection. Extracellularly recorded spikes were detected using the PTSD 
(Precise Timing Spike Detection) algorithm*^ Briefly, spike trains were built using 
three parameters: (1) a differential threshold set to 8 times the standard deviation of 
the baseline noise independently for each channel; (2) a peak lifetime period (set at 
2 ms); (3) a refractory period (set at 1 ms). The data presented in the text were not 
spike sorted. This choice was made according to the fact that, both in 2D and 3D 
networks, during a burst a global increase of the activity produces a fast sequence of 
spikes with different and overlapping shapes which make the sorting difficult and 
unreliable^^. Burst detection. Bursts were identified according to the method described 
in Pasquale et al.*". The algorithm is based on the computation of the logarithmic 
inter-spike interval histogram in which inter-burst activity (i.e., between bursts and/ 
or outside bursts) and intra-burst activity (i.e., within burst) for each recording 
channel can be easily identified, and then, a threshold for detecting spikes belonging 
to the same burst is automatically defined. Network burst detection. A similar 
procedure was followed for the detection of network bursts, looking for sequences of 
closely-spaced single- channels bursts. A network burst is identified if it involves at 
least the 80% of the network active channels*". Instantaneous Firing Rate (IFR). The 
IFR is computed by dividing the number of spikes which falls in a small sliding 
window by the window length. Such a window was implemented using a Gaussian 
kernel (100 ms width). The mean firing rate (MER) of the network was obtained by 
computing the firing rate of each channel averaged among all the active electrodes of 
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the MEA (i.e, being an active electrode defined by an activity > 0.1 spike/s). Peri or 
Post Stimulus Time Histogram (PSTH). The PSTH was calculated by considering 2 s 
time windows from the recordings that follow each stimulus. First, each time window 
was divided into 10 ms bins and then the number of spikes occurring in each time bin 
was counted. Finally, the histograms of each channel were normalized over the 
number of stimuli (50), the bin size and the maximum number of evoked spikes in a 
bin for each experiment (culture). Cross-correlation. The level of synchronization 
among multi-unit recordings was estimated by the cross-correlation function*^. 
Given two spike trains {e.g., X and Y) recorded from two electrodes of the MEA, we 
counted the number of events in the Y train within a time frame around the X event of 
±1 s, by using bins of 5 ms. The cross-correlogram coefficient, C(0), represents the 
area of the cross-correlogram in the central bin and it was evaluated to quantify the 
synchronization level among the recording channels. We defined an arbitrary 
threshold as <C(0)> + 2 std(C(0)) to detect the most synchronized electrodes. To 
quantify the correlation level the Cpea^ value was computed as the value of the cross- 
correlogram in an area around the maximum detected peak. C(0) and Cpeak values are 
normalized with respect to the auto-correlation value. Statistics. Unless otherwise 
stated, data were expressed as mean ± standard error of the mean. Statistical analysis 
was performed using MATLAB (The Mathworks, Natick, MA, USA). Since data do 
not follow a normal distribution (evaluated by the Kolmogorov-Smirnov normality 
test), we performed a non-parametric Mann- Whitney U-test. Significance levels were 
set atp < 0.001. 

MEA set-up 

The electrophysiological activity was recorded for 30 min by means of Micro- 
Electrode Arrays (MEAs) made up of 60 planar microelectrodes (TiN/SiN, 30 [xm 
electrode diameter, 200 |im spaced) arranged over an 8 X 8 square grid (except the 
four electrodes at the corners), supplied by Multi Channel Systems (MCS, Reutlingen, 
Germany). After 1200X amplification (MEA 1060, MCS), signals were sampled at 
50 kHz using the MCS data acquisition card. Recordings were performed outside the 
incubator at temperature of 37'^C; to prevent evaporation and changes of the pH 
medium, a constant slow flow of humidified gas (5% CO2, 20% O2, 75% N2) was 
inflated onto the MEA. Electrical stimuli were delivered by using a commercial 
general- purpose stimulus generator (STG 1008, MCS), which supplies current and 
voltage pulses to be applied to selected electrodes of the MEA (bottom stimulation), or 
to a borosilicate glass micropipette (estimated impedance of 0.9 MQ.). A tungsten coil 
(Concentric Bipolar Electrode, World Precision Instruments) was then inserted into 
the micropipette and connected to the stimulus generator (top stimulation). Stimuli 
were sent from the electrode at frequency of 0.2 Hzfor5 min with amplitude of 1.5 V 
(bottom stimulation) and 3 V (top stimulation) peak-to-peak. The stimulus pulse was 
biphasic (positive first) and lasted for 500 |j.s with a 50% duty cycle. 
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